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Abstract 

The state of the atmosphere is governed by the classical laws of fluid 
motion and exhibits correlations in various spatial and temporal scales. 
These correlations are crucial to understand the short and long term 
trends in climate. Cirrus clouds are important ingredients of the atmo- 
spheric boundary layer. To improve future parameterization of cirrus 
clouds in climate models, it is important to understand the cloud prop- 
erties and how they change within the cloud. We study correlations in 
the fluctuations of radar signals obtained at isodepths of winter and fall 
cirrus clouds. In particular we focus on three quantities: (i) the backscat- 
tering cross-section, (ii) the Doppler velocity and (iii) the Doppler spectral 
width. They correspond to the physical coefficients used in Navier Stokes 
equations to describe flows, i.e. bulk modulus, viscosity, and thermal con- 
ductivity. In all cases we find that power-law time correlations exist with a 
crossover between regimes at about 3 to 5 min. We also find that different 
type of correlations, including 1/f behavior, characterize the top and the 
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bottom layers and the bulk of the clouds. The underlying mechanisms for 
such correlations are suggested to originate in ice nucleation and crystal 
growth processes. 

1 Introduction 

Ice research spans a broad range of subjects, including studies of extremely 
complex phase diagrams [Mishima and Stanley, 1998], and phase transition 
modelling [Huckaby et al. 2000], hot topics as still seen recently [Putrino and 
Parrinello, 2002] in ab initio calculations, or in various properties, such as 
the dielectric constant [Gingle and Knasel, 1975] or acoustic emission features 
[Weiss and Grasso, 1997] during crack propagation. In clouds, ice crystals 
appear in a variety of forms and shapes, depending [Rogers, 1976; Heymsfield 
and Piatt, 1984] on the formation mechanism and the atmospheric conditions. 
In cirrus clouds, at temperatures lower than about —40° C ice crystals form and 
exist as mainly nonspherical particles [Petrenko and Whitworth, 1999]. 

In air transport, ice covering of airplane wings is also known to be very 
annoying. Moreover high altitude clouds, found above ca. 4 km including 
stratosphere cirrus [Tabazadeh et al, 1997] at > 20 km, play an important 
role in regulating the energy budget [Wallace and Hobbs, 1977] of the earth 
atmosphere system through their interaction with solar and terrestrial radiation 
[Stephens, 1990]. Their impact on the climate and weather systems at various 
time and space scales is not well understood. Recent studies suggests that the 
water vapor itself rather than the greenhouse gases could be the culprit in the 
global warming [Maurellis, 2001; Rosenfeld and Woodley, 2001]. Indeed due to 
the presence of water everywhere and the ability to change its state from, say, 
its frozen form at the poles to its liquid and vapor states in the atmosphere, the 
water is an effective transfer energy medium around the globe. No need to say 
that the description of ice crystal motions in clouds, and the cloud motion itself 
imply solving Navier-Stokes equations containing physical coefficients not well 
determined at this time. 

The radiative properties of the ice particles are important ingredients of 
the parameterization algorithms [Heymsfield and Piatt, 1984; Sundqvist et al., 
1989; Zhang and McFarlane, 1995; McFarquhar and Heymsfield, 1997; Zurovac- 
Jevtic, 1999] used in radiative transfer schemes of cirrus clouds in large-scale 
climate models. The question is still open on how cirrus radiative transfer 
properties might be predicted given their (actually) roughly known microphysi- 
cal characterization [Smith and Del Genio, 2001]. Therefore, it is essential that 
correlations between the structure of cirrus clouds and the radiative properties 
be better understood. An interesting study on composite cirrus morphology 
and correlations with temperature and large-scale vertical motion has been also 
conducted in e.g. [Mace, 1997]. 

Because of the vertical extent, ca. from about 4 to 14 km and higher, and the 
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layered structure of such clouds, one way of obtaining some information about 
their properties is using ground-based remote sensing instruments. (Satellites 
can also be used [Liou et al., 2002]). Experienced workers at airport radar 
sites have some feeling about ice cloud content. However a scientific inner 
cloud morphology investigation consists in searching for the statistical prop- 
erties of radio wave signals backscattered from the ice crystals present in the 
cloud. Backscattered signals received at the radar receiver antenna are known 
to depend on the ice mass content and on the particle size distribution [Atlas 
et al., 1995]. Ice particles can induce modifications in the radar signals which 
lead to fluctuations in (i) the backscattering cross section, (ii) the Dopplcr ve- 
locity and (iii) the Doppler spectral width. The motivation for our study is to 
identify properties of the inner structure of the cirrus clouds through a statis- 
tical analysis of the fluctuations of such measured signals. This is a first step 
toward including these statistics in parameterization schemes of cirrus clouds in 
large-scale climate models (see references quoted supra) or in short scale Navier- 
Stokes equation solutions. It is of interest to examine the time correlations in 
the backscattered signals not only on the boundary layers, i.e. at the top and 
bottom of the cloud, but also at several levels within the cloud, because of the 
extended vertical structure of cirrus clouds. Analysis of spatial correlations in 
vertical direction will be considered in a future study. 

Due to the nonlinear physics laws governing the phenomena in the atmo- 
sphere, data series of atmospheric quantities are usually non — stationary. 
Therefore a mere study based on the power spectrum only would be incom- 
plete, if not misleading. New techniques have been recently developed that 
can systematically eliminate trends in the data and thus reveal some intrin- 
sic dynamical properties such as fluctuation correlations that are very often 
masked by nonstationarities. There have been indeed a large number of stud- 
ies on long-range power-law correlations in time series [Kantz and Schreiber, 
1997; Brockwell and Davis, 1998; Malamud and Turcotte, 1999; Schreiber, 1999] 
in many research fields such as biology [Stanley et al, 1993; Hausdorff et al., 
1997; Ivanov et al, 1996, 1999; Mercik et al., 2000; Ashkenazy et al., 2001], 
finance [Mantegna and Stanley, 1995; Vandewalle and Ausloos, 1997; Ausloos 
and Ivanova, 1999] meteorology and climatology [Koscielny-Bunde et al, 1993, 
1998; Davis et al, 1996; Ivanova et al., 2000, 2002; Ausloos and Ivanova, 2001; 
Bunde et al, 2001; Peters et al, 2002 ]. Some of us have already examined 
long-range time correlations in nonstationary atmospheric signals of the liquid 
water path [Ivanova et al., 2000] and water vapor path [Ivanova et al., 2002] in 
stratus clouds, as well as in large-scale meteorological signals like the Southern 
Oscillation Index [Ausloos and Ivanova, 2001]. 

We pursue the analysis of millimeter wave cloud radar data [Mace et al., 
1997, 1998], hereby collected from observations during two consecutive winter 
days, i.e. Jan. 26 and 27, 1997 and one fall day, i.e. Sept. 26, 1997 at the 
Southern Great Plains site of the Atmospheric Radiation Measurements (ARM) 
program. Mace et al., [1998] have derived the microphysical properties of cir- 
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rus layer from surface based millimeter radar and infrared radiometer data. 
Mace et al. [1998] have developed techniques for retrieving cirrus cloud bound- 
aries, ice water content, crystal effective radius, and number concentration from 
reflectivities combined with atmospheric emitted radiance interferometer down- 
welling radiances and have applied the algorithm to almost two years of MMCR 
data (11/96-5/98) at the Atmospheric Radiation Measurement (ARM) Program 
Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site. Here 
we analyze the statistical properties of the three quantities: (i) the backscatter- 
ing cross section per unit volume, (ii) the Doppler velocity and (iii) the Dopplcr 
spectral width, obtain their time correlations, and give some physical interpre- 
tation. 

First, let Af(D,r) be the distribution function for the number of particles 
per unit volume; each particle is supposed to have a characteristic size D when 
located at r; let a(D) be the backscattering cross section of particles with such 
a characteristic size D. Then the total backscattering cross section (BCS) per 
unit volume [Clothiaux et al., 1995] is 

77(f) = [ a(D)J\f(D, r)dD. (1) 
Jd 

Leaving aside technical factors, the power received at the radar receiver antenna 
is proportional to 77(f) 

P(f) ~ ?W. (2) 

Due to the particle movement, frequency shifts w m , called Doppler shifts, appear 
in the scattered signal 

Wrn = " , (3) 

where m is the index of the particle that moves with its so called fall velocity 
v m and A c is the wavelength of the sounding signal [Clothiaux et al., 1995]. 

The radar is equipped with signal processing units, which strip off the carrier 
frequency and create signal voltages. By applying a Fourier transform over the 
raw voltage time series one can obtain the power density S m at each velocity v m 
[Clothiaux et al., 1995]. From Eq. (3) one can express the power density as a 
function of velocities rather than as a function of frequency. As such, the power 
density S m at each velocity v m is directly related to the total backscattering 
cross-section of all particles moving with a radial velocity v m with respect to 
the radar. Let S m be the normalized power density spectra at velocity v m , i.e. 

Sm = (4) 

where M corresponds to the number of particles in the resolution volume. 
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The radar signal processor units store the mean power-weighted radial ve- 
locity within the resolution volume as 

N 

< v >=^2 v m S m , (5) 

m— 1 

which is usually referred to as the Doppler velocity. 

The standard deviation of the power-weighted velocities about the mean is 
also stored 

AT 

°l = J2[v m ~ < v >} 2 S m (6) 

m— 1 

and is called the Doppler spectral width. These will be related to physical 
properties later on, in Section III. 

We study the signals of the backscattering cross section rj, the Doppler veloc- 
ity < v > and the Doppler spectral width <7% at different levels within the vertical 
structure of cirrus clouds. We apply the detrended fluctuation analysis (DFA) 
method [Peng et al, 1994] to characterize the correlations in these signals. The 
DFA method is suited to accurately quantifying power-law correlations in noisy 
nonstationary signals with polynomial trends [Peng et al, 1994; Vandewalle and 
Ausloos, 1998; Hu et al, 2001; Chen et al, 2002]. The advantage of the DFA 
method over conventional methods, such as the power spectrum analysis, is that 
it avoids the spurious detection of apparent long-range correlations that are an 
artifact of the nonstationarity (related to linear and higher order polynomial 
trends in the data). The DFA method has been tested at length by some of us 
on controlled time series that contain long-range correlations superposed on a 
(nonstationary) trend. Of note is a recent independent review of fractal fluctu- 
ation analysis methods which determined that DFA is one of the most robust 
methods [Taqqu et al, 1995]. 

Briefly, the DFA method involves the following steps: 

(i) The signal time series £(i), i — 1,2, ...,7V is integrated, to "mimic" a 
random walk 

2/(n) = £? =1 (£(i)-<£>), (7) 

where < £ >= Ej^f 

(ii) The integrated time series is divided into boxes of equal length, r. 
(in) In each box of length r, a least squares line is fit to the data (representing 

the trend in that box) . The y coordinate of the straight line segments is denoted 
by z(n). 

(iv) The integrated time series, y(n), is detrended by subtracting the local 
trend, z(n), in each box. 

(v) For a given box size r, the characteristic size of fluctuation F(t) for this 
integrated and detrended time series is calculated: 
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(vi) The above computation is repeated over all time scales (box sizes r) 
to provide a relationship between F(t) and the box size r (i.e. the size of the 
window of observation). 

A power law relation between the average root-mean-square fluctuation func- 
tion F(t) and the size of the observation window indicates the presence of so 
called scaling [Ma, 1976], i.e. when the fluctuation correlations can be charac- 
terized by a scaling exponent a, the self-similarity parameter, as in 



F(r)~T a . (9) 

For uncorrelated random walk fluctuations the F(t) function scales with an 
exponent a = 1.5 [Schroeder, 1991; Addison, 1997, Turcotte, 1997]. An exponent 
a > 1.5 indicates a persistent correlated behavior in the fluctuations, while an 
exponent a < 1.5 characterizes anticorrelations. In other words, if a < 1.5, the 
signal increments during the nonoverlapping successive time intervals of size r 
tend to be of opposite signs, so that the signal y(n) has a tendency to decrease in 
the next time step if it has had an increasing tendency in the previous one and 
vice versa. The feature is called antipersistency, or there are anticorrelations. 
If a > 1.5 the signal increments tend to have the same signs, so that y(n) tends 
to increase in the future if it has had an increasing tendency in the past, and 
conversely. The feature is called persistency. Physically, a persistent system is 
going to increase its statistical mean with time as if a positive feedback dom- 
inates the system. The antipersistency expresses a tendency of the values of 
increments to compensate for each over; this prevents the process from blowing 
up or going down very fast. Such a system tends to eliminate deviations, as 
if there is a negative feedback. An antipersistent time series visit, on average, 
the mean value more often than an ordinary one; the signal is said to present 
anticorrelations. Notice that the scaling exponent a can be related to the ex- 
ponent (3 of the power spectral density S(f ) <~ / _/3 of such a time series y(t), 
i.e. (5 = 2a — 1 [Heneghan and McDarby, 2000]. 

There are different ways of detrending the signal: one can fit the data in a 
window to a linear, quadratic, cubic [Vandewalle and Ausloos, 1998] or higher 
degree polynomial and then study the fluctuations of the signal. Discussions on 
different ways of detrending a signal, including a periodic and power-law trends, 
can be found in Hu et al. [2001]. 
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2 Data and Data Analysis 



We analyze data collected with a millimeter wave cloud radar that operates 
at 34.86 GHz ± 200 MHz. 1 The radar produce measurements at four modes 
[Clothiaux et at, 2000] that are used to obtain a best estimate to cloud signal 
following an interpolation procedure described in [Clothiaux et at, 2000]. The 
radar measures the backscattering cross-section of the emitted wave per unit 
volume of cloud particles. Assuming that the Rayleigh approximation is valid 
and ice crystals are spherical particles (we are aware of the approximation) this 
cross section is proportional to the sixth moment of the particle size and also 
depends on its ice water content [Mace, 1997]. The logarithm of the backscat- 
tering cross section is called the radar reflectivity and is the quantity typically 
used in the field. However, because we aim at studying the fluctuations, we 
prefer to use the raw backscattering cross section signal. 

Continuous measurements of (i) the backscattering cross section per unit 
volume n, (ii) the Dopplcr velocity < v > and (iii) the Doppler spectral width 
<7^ are recorded with 10 s temporal resolution 24 hours a day and with a 45 m 
spatial vertical resolution at 512 height levels. Radar reflectivity data measured 
with an accuracy of 0.5 dB are shown in Fig. 1. 

An order of magnitude value for the horizontal velocity at which such clouds 
move is 10 m/s. Therefore the volume examined at each data point can be 
considered to reflect the local ice content and velocity fluctuation in a fixed 
volume, as if the cloud is locally stable. The number of particles in the volume 
is assumed to be a constant during the measurements. If this is not the case, 
from a statistical physics point of view, one should work in a grand canonical 
ensemble formalism and subsequently introduce another ad hoc model than in 
[Clothiaux et at, 2000]. 

Studies along isotherms or isobars would be of great interest. This informa- 
tion being missing, we have decided to cut the cloud thickness into isodepths, 
identifying the cloud as being between its top and its bottom layer, i.e. at which 
the signal level is above -45 dB. We study the time series of the three sets of 
measurements at certain isodepths (relative to the thickness h of the cloud), i.e. 
top, 0.75 h, 0.5 h, 0.25 h, and bottom. These profiles are represented in Figs. 2 
and 3, from which the cloud height, thickness and their time evolution can be 
read. 

In this study in order not to introduce spurious time correlation effects, we 
avoid large data gaps within the cloud life time, e.g. the gaps observed between 
9:00 and 18:00 GMT on Jan. 26, 1997. Thus we consider that the time extent 
of this winter cloud at this location was between 18:00 GMT on Jan. 26 and 
06:00 GMT on Jan. 27, 1997, which amounts to 4315 data points. The fall 
cloud extends from 18:00 to 24:00 GMT on Sept. 26, 1997 which amounts to 
2048 data points. 

http: / /www.arm.gov/docs / sites/sgp/sgp.html 
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The time series of the backscattering cross section rj, Doppler velocity < v >, 
and Doppler spectral width cr 2 at the specified heights (or isodepths) within the 
winter cirrus cloud are plotted in Fig. 4(a-c). Results for each DFA-function 
F(t) are plotted in Fig. 5(a-e). The DFA- functions are grouped depending on 
the profile level, i.e. they refer to the top, 0.75 h, 0.5 h, 0.25 h, or the bottom 
of the cloud. For the fall cirrus cloud (observed at the same experimental site) , 
the data series of the three quantities of interest have been also analyzed along 
the same profile levels as the winter cloud (see Fig. 6). The results from the 
DFA analysis at the different height levels are displayed in Fig.7(a-e). 

The scaling analysis has been used performing statistical tests to find the 
best slopes for the fits, and the corresponding scaling regimes. Each scaling 
behavior is characterized by slightly different scaling exponents, cti and ct2 at 
low and large time lag r respectively; we find a crossover at t x ~ 3 — 5 min 
for all data series and physical quantities. The values of the a-exponents for 
the two scaling regimes and for the three studied quantities are summarized 
in Table 1 and Table 2 for the winter cloud and the fall cloud. Their values, 
whence the clouds, can be compared with each other, from signal type to signal 
type, and on both sides of the crossover time lag for a given signal. The error of 
the estimate of each a value is thus the error of the slope and is defined as the 
square root of residual sum of squares divided by the number of points in the fit 
minus two and by the sum of squares of the deviation with respect to the mean 
of the time lags [Bowman and Robinson, 1987]. The correlation coefficient for 
all estimates is in the range between 0.992 and 0.998. 

We find that the a\ values have a maximum in the center of the cloud for 
both cloud cases. We note an increase of a\ from bottom to top for the Doppler 
velocity of the winter cloud, while the value of a.\ is rather constant for the fall 
cloud. The a.\ values remain constant for the 77, < v > and cr 2 data and for 
both cloud cases, and close to 1.5, except for < v > at the cloud bottom in the 
Jan. 26 and 27, 1997 case. 2 

The «2 values are systematically close to or below 1.0, for < v > and <7 2 . 
The backscattering cross section 77 values fall around 1.0. The 012 values decrease 
from bottom to top for < v > and cr 2 , but present a maximum in the fall cloud 
center for the r\ data. The U2 values reach a smooth maximum in the center of 
the cloud for both < v > and cr 2 data. The a.2 values reach a maximum in the 
center of the cloud for the 77, < v > and cr 2 data and so in both cloud cases, 
though all these values are quite close to 1.0. 

Thus, scaling exponents at the top of the cloud seem to have usually a lower 
value as compared to that in the bulk of the cloud, - except for < v >. This could 
be due to a physical effect or arise from the instrumental noise contribution: at 

2 The short-range correlations may be influenced by the interpolation procedure that pro- 
duces the best estimate to cloud signal from the four modes measurements of the millimeter 
wave radar. The current interpolating procedure is complicated [Clothiaux et ai, 2000] and 
mimics the signal from the four modes in a systematic but non-uniform way. 
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such heights the real signal may be low as compared to the instrumental noise. 3 
To check the robustness of the results it is essential to perform a test on 
surrogate data [Schreiber, 1999; Schreiber and Schmitz, 2000]: we randomly 
shuffled the amplitudes of (i) the backscattering cross section per unit volume 
r\, (ii) the Doppler velocity < v > and (iii) the Doppler spectral width er^ signals 
for winter and fall clouds. It is known [Viswanathan et al, 2002] that the fat 
tailed distributions are thought to be caused by long-range volatility correla- 
tions. Destroying all correlations by shuffling the order of the fluctuations, is 
known to cause the fat tails almost to vanish. We have found that the long- 
range correlations do vanish in surrogate data as seen from the DFA-functions 
plotted in Figs. 8 and 9 and Tables 3 and 4. This, together with an excellent 
systematics in the root mean square linear fit correlation coefficients, leads to 
put substance to the results found here above. 

3 Discussion 

Four main points seem to need some emphasis and bring newly interesting 
results: 

(I) We find that DFA-functions F(t) for all radar backscattering time series 
exhibit a crossover at about 4 min. We observe a different scaling behavior for 
the backscattering cross section n, the Doppler velocity < v > and the Doppler 
spectral width a% between the top of the cloud from the bulk, and from the 
bottom, especially for large time scales. Short time range (< 4 min) scaling 
of n is completely uncorrelated with ot\ « 1.5. Long-range correlations of n, 
< v > and o~y are quite similar (within error bars) for the bulk of the cloud with 
a 2 ~ 0.90 for < v > and and a 2 ~ 1.00 for n. Such a-exponents (or (3=1) 
characterize a special physical phenomenon of interest [Schroeder, 1991], that is 
1/f noise [Shlesinger, 1987; Weismann, 1988; Shlesinger and West, 1988; West 
and Shlesinger, 1989]. 

(II) Cirrus cloud temperature measurements with radiosondes show that 
the temperature within the cloud field usually varies smoothly between about 
—43° C and —33° C. Spontaneous freezing is thought to occur at temperatures 
colder than about —40° C. Assuming therefore that kinetically driven homo- 
geneous nucleation is the primary formation mechanism at these temperatures 
[Kiang et. al, 1971; Sassen and Dodd, 1988; Tabazadeh et al., 1997a, 1997b] 
supercooled liquid water drops are thus seen to nucleate (and die rapidly) in 
water-satured updrafts. 

(III) However beside the a — 1.5 value, i.e. indicating no correlation at all in 
the signal fluctuations at short time ranges, thus suggesting a birth and death 
(nucleation) process, the scaling difference at the crossover, to a 1/f process, 
suggests propagation effects in the cirrus cloud at long time scales, similar to 
growth and percolation features. The Doppler velocity measurements of the 

3 http : / /www. arm. gov /docs /instruments /static/ 'mmcr.html 
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cirrus ice particles shown in Fig. 2b suggest updrafts in fact, ... therefore 
updrafts serve as mechanisms for the aggregation of ice particles at long time 
ranges. The data analysis is consistent with (and somewhat proves) such an 
intuitive physical mechanism. 

Further physical insight is obtained if we translate the n, < v > and al 
data from scattering information into usual fluid mechanics parameters. Even 
though the data is analyzed through a model with constant number of particles 
with a given (spherical) shape, and other standard approximations, neverthe- 
less the backscattering cross section is a measure of the scattering intensity, thus 
clearly a measure of the ice content, in principle including particle shape, size 
and distribution. Therefore the ?7-DFA function is the direct signature of the 
density-density correlation function, i.e. the zero-wave vector structure factor. 
In other words, it is a measure of the cloud (fluid) compressibility or bulk mod- 
ulus [Chaikin and Lubensky, 1995]. The < v > data leads to information about 
the particle velocity in the cloud, while the < v >-DFA function relates to the 
velocity-velocity correlation function. Even assuming a constant density, or a 
constant and identical mass for each particles, (... we are aware that this is an 
approximation), < v >-DFA measures the correlations in particle momenta and 
thus relates to the viscosity [Huang, 1967] of the inner cloud structure. On the 
contrary, if the particle density is not homogeneous the < v >-DFA function 
still has a physical meaning, i.e. it is a signature for the diffusivity or diffusion 
coefficient [Chaikin and Lubensky, 1995]. Finally the data is a measure of the 
dissipated power due to inelastic scattering, through the velocity square [Huang, 
1967], and is thus related to the (cloud) temperature. The cr^-DFA function is 
thus a measure of the temperature-temperature correlation function, whence the 
thermal conductivity [Huang, 1967]. 

In principle, the equations of motion for the above three physical quantities, 
i.e. ice density, ice crystal velocity and temperature in such cirrus clouds can 
be written as Navier-Stokes equations, - in terms of these three basic physi- 
cal coefficients. They constitute an essential knowledge to be acquired before 
meaningful physical modelling. Moreover, the diffusion coefficient and the vis- 
cosity coefficient can serve as inputs in the corresponding microscopic evolution 
equations, i.e. the Langevin equations, with appropriate noise and drift terms. 

In conclusion, the time dependence characteristics of physical quantities in 
cirrus clouds have been obtained. Short time range correlations of radar back 
scattering cross section are uncorrelated, i.e. are Brownian-like; at time lags 
longer than about 3 to 5 min correlations are of the 1// noise type. Cer- 
tain asymmetry is seen to exist from cloud top to bottom in respective scaling 
properties. This can be understood to be due to different mechanisms on the 
boundaries and in the bulk of the cloud, depending on the time scales. The 
data analysis leads to a consistent physical picture on the inner noncquilibrium 
structure, i.e. ice crystal nucleation and growth. Finally, the physical descrip- 
tion of such cirrus clouds, is shown to be related to the bulk modulus, viscosity 
and thermal conductivity. 
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Figure Captions 

Figure 1 - Radar reflectivity observations on (a) Jan. 26, (b) 27, 1997 
and (c) on Sept. 26, 1997 at the Southern Great Plains site of Atmospheric 
Radiation Measurements program. 

Figure 2 - Cloud profiles at top and bottom and at 0.75, 0.5 and 0.25 of the 
thickness h of the winter cirrus cloud (Fig. la, b) along which the backscattering 
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cross section, Doppler velocity and Doppler spectral width are studied. 

Figure 3 - Cloud profiles at top and bottom and at 0.75, 0.5 and 0.25 of 
the thickness h of the fall cirrus cloud along which the backscattering cross 
section, Doppler velocity and Doppler spectral width are studied. Measure- 
ments are obtained at the Southern Great Plains site of Atmospheric Radiation 
Measurements program on Sept. 26, 1997 (Fig. lc). 

Figure 4 - (a) Backscattering cross section per unit volume at different 
relative heights to the thickness h of the cirrus cloud as measured on Jan. 
26 and 27, 1997 at the Southern Great Plains site of Atmospheric Radiation 
Measurements program (data in Fig. la,b). (b) Doppler velocity (m/s). The 
positive spike between 20 and 22 h of the "bottom" data series is of order of 
11 m/s. (not shown) at heights in the cloud shown in Fig. 2.; (c) Doppler spectral 
width (m/s) at heights in the cloud shown in Fig. 2. 

Figure 5 - DFA- functions for backscattering cross section (circles) , Doppler 
velocity (triangles) and Doppler spectral width (squares) at the (a) top, (b) 0.75, 
(c) 0.50, (d) 0.25 relative to the thickness h of the cloud, and (e) bottom of the 
winter cloud shown in Fig. 2. DFA-functions are displaced for readability. The 
a- values are listed in Table 1. 

Figure 6 - (a) Backscattering cross section per unit volume at different 
relative heights to the thickness h of the cirrus cloud as measured on Sept. 26, 
1997 at the Southern Great Plains site of Atmospheric Radiation Measurements 
program (Fig. lc). (b) Doppler velocity (m/s) at heights in the cloud shown in 
Fig. 3; (c) Doppler spectral width (m/s) at different heights in the cloud shown 
in Fig. 3. 

Figure 7 - DFA-functions for backscattering cross section (circles), Doppler 
velocity (triangles) and Doppler spectral width (squares) at the (a) top, (b) 0.75, 
(c) 0.50, (d) 0.25 relative to the thickness h of the cloud, and (e) bottom of the 
cloud shown in Fig. 3. DFA-functions are displaced for readability. The a- 
values are listed in Table 2. 

Figure 8 - DFA-functions for shuffled data of: backscattering cross section 
(circles), Doppler velocity (triangles) and Doppler spectral width (squares) at 
the (a) top, (b) 0.75, (c) 0.50, (d) 0.25 relative to the thickness h of the cloud, and 
(e) bottom of the winter cloud shown in Fig. 2. DFA-functions are displaced 
for readability. The a- values are listed in Table 3. 

Figure 9 - DFA-functions for the shuffled data of: backscattering cross sec- 
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tion (circles) , Doppler velocity (triangles) and Doppler spectral width (squares) 
at the (a) top, (b) 0.75, (c) 0.50, (d) 0.25 relative to the thickness h of the cloud, 
and (e) bottom of the fall cloud shown in Fig. 3. DFA- functions are displaced 
for readability. The a- values are listed in Table 4. 
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Table 1: Values of the a\ and a2-exponents from the DFA analysis of the 
backscattered cross section n, Doppler velocity < v > and Doppler spectral 
width ai data shown in Fig. 4(a-c); winter cloud. The crossover a\ to a 2 is 
at about 3 min. The error of the estimate of each a value is thus the error of 
the slope and is defined as the square root of residual sum of squares divided 
by the number of points in the fit minus two and by the sum of squares of the 
deviation with respect to the mean of the time lags [Bowman and Robinson, 
1987]. 



position 


V 


< v > 






ai 


a 2 


Oil 


a.2 


ai 


a.2 


top 
0.75 h 
0.50 h 
0.25 h 
bottom 


1.36 ±0.03 
1.41 ±0.04 
1.62 ±0.10 
1.43 ±0.04 
1.46 ±0.09 


0.89 ±0.01 
1.11 ±0.01 
0.94 ±0.01 
1.15 ±0.02 
1.27 ±0.03 


1.43 ±0.07 
1.38 ±0.04 
1.42 ±0.04 
1.29 ±0.03 
1.09 ±0.05 


0.76 ±0.01 
0.88 ±0.01 
0.96 ±0.01 
0.93 ±0.01 
0.70 ±0.02 


1.24 ±0.04 
1.36 ±0.03 
1.35 ±0.04 
1.31 ±0.04 
1.38 ±0.05 


0.73 ±0.01 
0.90 ±0.02 
0.92 ±0.02 
0.90 ±0.02 
1.02 ±0.01 



Table 2: Values of the a\ and (^-exponents from the DFA analysis of backseat - 
tering cross section n, Doppler velocity < v > and Doppler spectral width 0% 
data shown in Fig. 6(a-c); fall cloud. A crossover is observed at about 5 min. 



position 


V 


< v > 


< 






a 2 


Oil 




Oil 


a 2 


top 
0.75 h 
0.50 h 
0.25 h 
bottom 


1.32 ±0.05 
1.63 ±0.04 
1.50 ±0.04 
1.69 ±0.04 
1.36 ±0.07 


0.66 ±0.04 
0.85 ±0.05 
1.13 ±0.04 
1.23 ±0.04 
0.70 ±0.03 


1.34 ±0.06 
1.14 ±0.05 
1.14 ±0.03 
1.31 ±0.04 
1.22 ±0.04 


0.66 ±0.01 
0.86 ±0.03 
0.85 ±0.03 
0.93 ±0.03 
0.68 ±0.03 


1.32 ±0.03 
1.37 ±0.03 
1.37 ±0.03 
1.27 ±0.04 
1.22 ±0.03 


0.69 ±0.03 
1.07 ±0.05 
1.19 ±0.05 
1.09 ±0.03 
0.88 ±0.02 
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Table 3: Values of the a-exponents- from the DFA analysis of the shuffled data 
of backscattered cross section n, Doppler velocity < v > and Doppler spectral 
width &y data shown in Fig. 4(a-c); winter cloud. The correlation coefficient 
for all cases is around 0.997. 



position 


V 


< v > 


< 




a 


a 


a 


top 


0.47 ±0.006 


0.49 ±0.006 


0.51 ±0.003 


0.75 h 


0.53 ±0.004 


0.53 ±0.004 


0.54 ±0.004 


0.50 h 


0.52 ±0.004 


0.50 ±0.005 


0.52 ±0.005 


0.25 h 


0.56 ±0.006 


0.53 ±0.005 


0.48 ±0.006 


bottom 


0.49 ±0.006 


0.50 ±0.004 


0.49 ±0.005 



Table 4: Values of the a-exponcnt from the DFA analysis of shuffled data of 
backscattering cross section n, Doppler velocity < v > and Doppler spectral 
width Cy data shown in Fig. 6(a-c); fall cloud. The correlation coefficient for 
all cases is around 0.994. 



position 


V 


< v > 






a 


a 


a 


top 
0.75 h 
0.50 h 
0.25 h 
bottom 


0.55 ±0.006 
0.48 ± 0.006 
0.50 ±0.007 
0.49 ± 0.006 
0.51 ±0.010 


0.48 ±0.006 
0.52 ±0.007 
0.51 ±0.006 
0.52 ±0.005 
0.51 ±0.005 


0.49 ±0.005 
0.52 ±0.006 
0.46 ±0.006 
0.52 ±0.007 
0.53 ±0.007 
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Figure 1: Radar reflectivity observations on (a) Jan. 26, (b) 27, 1997 and (c) 
on Sept. 26, 1997 at the Southern Great Plains site of Atmospheric Radiation 
Measurements program. 
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Figure 2: Cloud profiles at top and bottom and at 0.75, 0.5 and 0.25 of the 
thickness h of the winter cirrus cloud (Fig. la, b) along which the backscattering 
cross section, Doppler velocity and Doppler spectral width are studied. 
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Figure 3: Cloud profiles at top and bottom and at 0.75, 0.5 and 0.25 of the 
thickness h of the fall cirrus cloud along which the backscattering cross sec- 
tion, Doppler velocity and Doppler spectral width are studied. Measurements 
are obtained at the Southern Great Plains site of Atmospheric Radiation Mea- 
surements program on Sept. 26, 1997 (Fig. lc). 
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Figure 4: (a) Backscattering cross section per unit volume at different relative 
heights to the thickness h of the cirrus cloud as measured on Jan. 26 and 27, 
1997 at the Southern Great Plains site of Atmospheric Radiation Measurements 
program (data in Fig. la,b). 
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(b) Doppler velocity (m/s) 
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Figure 4: (b) Doppler velocity (m/s). The positive spike between 20 and 22 h 
of the "bottom" data series is of order of 11 m/s (not shown) at heights in the 
cloud shown in Fig. 2 
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(c) Doppler spectra width (m/s) 
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Figure 4: (c) Doppler spectral width (m/s) at heights in the cloud shown in Fig. 
2 
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Figure 5: DFA- functions for backscattering cross section (circles), Doppler ve- 
locity (triangles) and Doppler spectral width (squares) at the (a) top, (b) 0.75, 
(c) 0.50, (d) 0.25 relative to the thickness h of the cloud, and (e) bottom of the 
winter cloud shown in Fig. 2. DFA-functions are displaced for readability. The 
a- values are listed in Table 1. 
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Figure 6: (a) Backscattering cross section per unit volume at different relative 
heights to the thickness h of the cirrus cloud as measured on Sept. 26, 1997 at the 
Southern Great Plains site of Atmospheric Radiation Measurements program 
(Fig. lc). 



28 



(b) Doppler velocity (m/s) 




18 19 20 21 22 23 24 



time (hours) 

Figure 6: (b) Doppler velocity (m/s) at heights in the cloud shown in Fig. 3. 
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(c) Doppler spectra width (m/s) 
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Figure 6: (c) Doppler spectral width (m/s) at heights in the cloud shown in Fig. 
3. 
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Figure 7: DFA- functions for backscattering cross section (circles), Doppler ve- 
locity (triangles) and Doppler spectral width (squares) at the (a) top, (b) 0.75, 
(c) 0.50, (d) 0.25 relative to the thickness h of the cloud, and (e) bottom of 
the cloud shown in Fig. 3. DFA-functions are displaced for readability. The 
a- values arc listed in Table 2. 
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Figure 8: DFA-functions for shuffled data of: backscattering cross section (cir- 
cles), Doppler velocity (triangles) and Doppler spectral width (squares) at the 
(a) top, (b) 0.75, (c) 0.50, (d) 0.25 relative to the thickness h of the cloud, and 
(e) bottom of the winter cloud shown in Fig. 2. DFA-functions are displaced 
for readability. The a- values are listed in Table 3. 
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Figure 9: DFA-functions for the shuffled data of: backscattering cross section 
(circles), Doppler velocity (triangles) and Doppler spectral width (squares) at 
the (a) top, (b) 0.75, (c) 0.50, (d) 0.25 relative to the thickness h of the cloud, 
and (e) bottom of the fall cloud shown in Fig. 3. DFA-functions are displaced 
for readability. The a- values are listed in Table 4. 
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